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Abstract 

The calculations in Thomas-Fermi approximation show that in a 
gravitational field each cell of ultra dense matter inside celestial bod- 
ies obtains a very small positive electric charge. A celestial body is 
electrically neutral as a whole, because the negative electric charge 
exists on its surface. On the order of magnitude the positive volume 
charge is very small (10 _18 e only). But it is sufficient to explain the 
occurrence of magnetic fields of celestial bodies and the existence of a 
discrete spectrum of steady-state values of masses of stars and pulsars. 
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We cannot measure magnetic fields of the majority of stars, which are distant 
far from us. Therefore, the existence of magnetic fields for the majority of 
stars can be considered only hypothetically. However, the magnetic field of 
the Sun is known over than a hundred years, and in the last decades the 
astronomers managed to measure magnetic fields for a number of stars (so- 
called Ap-stais) []I| and some pulsars 0. It is interesting to construct a model 
describing the generation of magnetic fields by stars and to compare it with 
the data of the astronomers. The mechanism examined below is based on 
the gravity- induced electric polarization of matter. It is capable to explain 
also the generation of magnetic fields by planets ||, however, in the case of 
stars, this mechanism works in the purest manner. 



The action of gravity on metals has often been a topic of discussion before 
[§]-@. The basic result of these researches is reduced to the statement that 
gravity induces inside a metal an electric field with an intensity 

m iQ , x 

E~-f-, (1) 

where rrii is the mass of an ion, 
g is gravity acceleration, 
e is the electron charge. 

This field is so small that it is not possible to measure it experimentally. 
It is a direct consequence of the presence of an ion lattice in a metal. This 
lattice is deformed by gravity and then the electron gas adapts its density to 
this deformation. The resulting field becomes very small. 

Under superhigh pressure, all substances transform into ultradense matter 
usually named nuclear-electron plasma [[31]]. It occurs when external pressure 



enhances the density of matter several times [ID, TT|. Such values of pressure 
exist inside celestial bodies. 

In nuclear-electron plasma the electrons form the degenerated Fermi gas. 
At the same time, the positively charged ions form inside plasma a dense 



packing lattice [12|,[T3[. As usually accepted, this lattice may be replaced 
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by a lattice of spherical cells of the same volume. The radius r s of such a 
spherical cell in plasma of the mass density 7 is given by 



£?=m~'=^ (2) 

3 \rriiJ n 

where Z is the charge of the nucleus, rrii = Am v is the mass of the nucleus, 
A is the atomic number of the nucleus, m p is the mass of a proton, and n is 
the electron number density 

The equilibrium condition in matter is described by the constancy of 
its electrochemical potential flOf . In plasma, the direct interaction between 
nuclei is absent, therefore the equilibrium in a nuclear subsystem of plasma 
(at T = 0) looks like 

fii = rriiip + Zap = const. (4) 

Here <p is the potential of an electric field and ip is the potential of a 
gravitational field. 

The direct action of gravitation on electrons can be neglected. Therefore, 
the equilibrium condition in the electron gas is 

He = 7^- - (e - 5q)(f = const, (5) 

where m e is the mass of an electron and pp is the Fermi momentum. 

By introducing the charge Sq, we take into account that the charge of the 
electron cloud inside a cell can differ from e. A small number of electrons can 
stay on the surface of a plasma body where the electric potential is absent. 
It results that the charge in a cell, subjected to the action of the electric 
potential, is effectively decreased on a small value Sq. 

The electric polarization in plasma is a result of changing in density of 
both nuclear and electron gas subsystems. The electrostatic potential of the 
arising field is determined by the Gauss' law 
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where the position of nuclei is described by the function 5(r). 
According to the Thomas - Fermi method, n is approximated by 



n 



8n 3 



(7) 



With this substitution, Eq.(|B|) is converted into a nonlinear differential 
equation for tp, which for r > is given by 



1 d ( 2 d . . 
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It can be simplified by introducing the following variables JLO 

,u 



H e + (e — 5q)(p = Ze 



(8) 



(9) 



and r = ax, 
where 

with an = — 2 = 



Bohr radius. 



With the account of Eq. ( 
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Then Eq.(H) gives 



d 2 u it 3//2 
dx 2 
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In terms of u and x, the electron density within a cell is given by [|Kj 
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Under the influence of gravity the charge of the electron gas in a cell 
becomes equal to 



Q e = 4ire / n{r)r dr = — — 
Jo 3/r 
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Using Eq.(O), we obtain 
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At r — > the electric potential is due to the nucleus alone tp(r) — 
means that u(0) — > 1 and each cell of plasma obtains a small charge 
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(14) 
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For a cell placed in a point R inside a star 



8q — Zeri 
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Considering that gravity acceleration g = — ^ and the electric field in- 



tensity E 



dip 
dR 



dr s 
dR 



2 rI fg + 8qE 
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This equation has the following solution 



dr s 

— = 
dR 



(17) 



(18) 



and 



^g + SqE 



0. 



(19) 



In plasma, the equilibrium value of the electric field on nuclei according 
to Eq.([|) is determined by Eq.(|I]) as well as in a metal. But there is one more 
additional effect in plasma. Simultaneously with the supporting of nuclei in 
equilibrium, each cell obtains an extremely small positive electric charge. 

As divg = —inGnmi and divE = 4iin5q, the gravity-induced electric 
charge in a cell 



Sq 



~Z ~ 



KT 18 e, 



(20) 



where G is the gravity constant. 
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However, because the sizes of bodies may be very large, the electric field 
intensity may be very large as well 



■a=± (2D 



In accordance with Eqs.(I5|,IPD, the action of gravity on matter is com- 
pensated by the electric force and the gradient of pressure is absent. 

Thus, a celestial body is electrically neutral as a whole, because the pos- 
itive volume charge is concentrated inside the charged core and the negative 
electric charge exists on its surface and so one can infer gravity-induced elec- 
tric polarization of a body. 



At the surface of the core, the electric field intensity reduces to zero. The 
jump in electric field intensity is accompanied at the surface of the core 
by the pressure jump Ap(R N ). It leads to the redistribution of the matter 
density inside a star. In a celestial body consisting of matter with an atomic 
structure, density and pressure grow monotonously with depth. In a celestial 
body consisting of electron-nuclear plasma, the pressure gradient inside the 
polarized core is absent and the matter density is constant. Pressure affecting 
the matter inside this body is equal to the pressure jump on the surface of 
the core 

P = A P (R N ) = = 2 -l G ^R%, (22) 

07T y 
where is the radius of the core. 

One can say that this pressure jump is due to the existence of the polar- 
ization jump or, which is the same, the existence of the bound surface charge 
formed by an electron pushed out from the core and making the total charge 
of the celestial body equal to zero. 

Because the electron subsystem of plasma inside a star is the relativistic 



Fermi gas, we can write its equation of state [10] 



P ~ 4 m// 3 /3 4 /3 [ ] 
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where f3 ■ m p is the mass of the matter related to one electron of the Fermi 
gas system, and 

m p is the proton mass. 

Because of the electroneutrality, one proton should be related to elec- 
tron of the Fermi gas of plasma. The existence of one neutron per proton is 
characteristic for a substance consisting of light nuclei. The quantity of neu- 
trons grows approximately to 1.8 per proton for the heavy nuclei substance. 
Therefore, it is necessary to expect that inside stars 

2 < (3 < 2.8. (24) 
As pressure inside a star is known (Eq.(|22|)), from Eq.(E3l) it is possible 



to determine a steady-state value of mass of a star 

M* = (Al^. (25) 
This mass is expressed by dimensionless constants only 

& = L54 ' 1038 (26) 

( = (1.5V) 1 / 2 ~ 5, 

and the slowly varying parameter j3 (Eq.(|24D). 

The masses of stars can be measured with a considerable accuracy, if 
these stars compose a binary system. There are almost 200 double stars 



whose masses are known with the required accuracy [|15|]. Among these stars 
there are giants, white dwarfs, and stars of the main sequence. Their averaged 
mass is described by the equality 

(M*) = (1.36 ± 0.05) M Q , (27) 
where M is the mass of the Sun. 

The center of this distribution (Fig.l) corresponds to Eq.(^5|) at (3 ~ 2.6. 

It is interesting to note that the "biography" of such a star appears much 
poorer than in the Chandrasecar model. 

Temperature does not influence the parameters of relativistic plasma. 
Therefore, a star with a mass close to the steady-state value (Eq.fl25|)) is in a 
stable equilibrium not depending on temperature. It should not collapse with 
a temperature decreasing. The instability of a star can arise with burning 
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Figure 1: Mass distributions of stars and pulsars from the binary systems 
The curve shows Eq.@. 
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out of light nuclei - deuterium and helium - and with a related increasing 
of j3. This growth leads to the reduction of a steady-state value of mass 
(Eq.(P5D) and, probably, to the distraction of stars with greater masses. 



As the density of matter inside a relativistic star is constant, it is possible 
to assume that it equals the mean density of the Sun and to estimate a star 
radius 

M, \V3 



R 

3 

where 7 Q is the mean density of the Sun. 



4?r / 

7©/ 



) , (28) 



It allows one to calculate the momentum of a star as the momentum of a 
sphere with a constant density 

I = -M±R 2 (29) 

and at a known frequency of rotation Q to calculate its angular momen- 
tum 

L = -M+QR 2 . (30) 

In this model the magnetic moment of a star is created by the rotation 
of a star as a whole. Thus, it is composed of two parts. One is the magnetic 
moment of the layer of electrons placed on the external surface of a star 

^ = -l{f P R 3 )nR 2 . (3i) 

The second component of the magnetic moment is created by the posi- 
tively charged core 

^ = 1{y pR3 ) qr2 - (32) 

The summary moment is 
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Figure 2: The observed values of the magnetic moments of celestial bodies 
vs. their angular momenta. On the ordinate, the logarithm of the magnetic 
moment over Gs ■ cm 3 is plotted; on the abscissa the logarithm of the angular 
momentum over erg-s is shown. The solid line illustrates Eq.(^). The dash- 
dotted line is the fitting of the observed values. 
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Ato = — (ypi^ntf. (33) 

It is remarkable that the gyromagnetic relation of a star, i.e., the relation 
of its magnetic moment to the angular momentum, is expressed through 
world constants only 

« = T = # (34 » 

The measurements permit us to define the frequency of rotation and mag- 
netic fields for a number of stars []]]. It appears enough to check up the 
considered theory, since masses of stars and their momenta are determined 
inside the theory (Eq.(|27D and (Eq.(|30|))). The magnetic moments as func- 
tions of their angular momenta for all celestial objects (for which they are 
known today) are shown in Fig.2. The data for planets are taken from [|TJ 



the data for stars are taken from [|l], and for pulsars - from p. As it can 
be seen from this figure with the logarithmic accuracy, all celestial bodies - 
stars, planets, and pulsars - really have the gyromagnetic ratio close to the 
universal value (Eq. (|34|)). Only the data for the Moon fall out, because its 
size is too small to create an electrically polarized core. 



Apparently, the considered theory is quite true for pulsars which consist, as 
it is supposed, from the neutron substance with an addition of electrons and 



protons [10]. As this substance is a relativistic one, there is a fair definition 



of a steady-state value of mass Eq.(|2B). The astronomers measured masses 



of 16 radio-pulsars and 7 x-ray pulsars included in a double system [[17 
According to this data, the distribution of masses of pulsars is 

(M ptljMr ) = (1.38±O.O3)M0. (35) 

The center of this distribution corresponds to Eq . (|25|) at (3 ~ 2.6. 

The gyromagnetic relations are measured for three pulsars only M . These 
values are in a quite satisfactory agreement with Eq. (|34]) (Fig.2). For the 
majority of pulsars |18|], there are estimations of magnetic fields obtained 



using a number of model assumptions Bfl. It is impossible to consider these 
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Figure 3: The estimated values of the magnetic moments of pulsars [|18j vs. 
their angular momenta. Solid line is Eq . fl34[) . Axes are like in Fig. 2. 
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data as the data of measurements, but nevertheless they also are in some 
agreement with Eq . (|34D , (Fig.3). 

For planets the situation is more difficult. First, inside planets the sub- 
stance forms not relativistic electron-nuclear plasma, but nonrelativistic electron- 
ion plasma. It has different equation of state fl(| leading to a more complex 
expression for the stable mass of a planet core than the expression of Eq.([25|) 
for stars. Second, a noncharged layer at the surface of the core can take a 
significant part of a planet's volume and it is impossible to neglect a role 
of this stratum. However, it can be seen from Fig. 2 that the gyromagnetic 
relations of planets are also in the quite satisfactory agreement with Eg.(|34|). 
The detailed calculation for the Earth || gives for the magnetic moment 
4 • 10 25 Oe • cm 3 , which is almost exactly twice smaller than the measured 
value 8.05 • 10 25 Oe ■ cm 3 . Thus, it is possible to assume that the basic com- 
ponent of the magnetic moment of planets is induced by the same mechanism 
which is working in stars. 



13 



References 

Bona E.F. and Landstreet J.D. - The Astrophysical Journ, Suppl., 1980, 
v.42, 421-445. 

Beskin V.S.,Gurevich, A.V., Istomin Ya.N. - Physics of the Pulsar Mag- 
netosphere, Cambridge University Press, 1993. 

Vasiliev B.V. - Nuovo Cimento B,1999,v.ll4,pp.291-300. 

Stuff L.I. and Barnhill M.V. - Phys.Rev.,1968,v.l51,pp.l067-1071. 

Dressier A.I. a.o. -Phys.Rev.,1968,v.l68,pp.737-743. 

Riegel T.J. - Phys. Rev.B,1970,v. 2,pp.825-828. 

Kumar N. and Naddini R. - Phys. Rev.D.,1973,v.7,pp.l067-1071. 

Leung M.C. et al. - Canad.Journ. of Phys.,1971,v.49,pp.2754-2767. 

Leung M.C. - Nuovo Cimento,1972,v.76,pp.825-929. 

Landau L.D. and Lifshits E.M. - Statistical Physics, 1980, vol.1, 3rd edi- 
tion, Oxford: Pergamon. 

Vasiliev B.V. and Luboshits V.L. - Physics-Uspekhi,1994,v.37,pp.345- 
351. 

Kirzhnitz D.A. - JETP, 1960, v.38, pp.503-508. 

Abrikosov A.A. - JETP, 1960, v.39, pp.1797-1805. 

Leung Y.C. - Physics of Dense Matter, 1984, Science 

Heintz W.D. - Double stars, 1978, Geoph. and Astroph. monographs, 
vol.15, D.Reidel Publ.Comp. 

Sirag S.-P. - Nature,1979,v.275,pp.535-538. 

Thorsett S.E. and Chakrabarty D. - E-preprint: |astro-ph/9803260| , 1998, 
35pp. 

[18] Taylor J.H., Manchester R.N., Lyne A.G., Camilo F., Catalog of 706 
pulsars, 1995, pulsar.prinston.edu 



14 



